Stabilization of colloidal suspensions by means of highly-charged nanoparticles 
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We employ a novel Monte Carlo simulation scheme to elucidate the stabilization of neutral col- 
loidal microspheres by means of highly-charged nanoparticles [V. Tohver et al., Proc. Natl. Acad. 
Sci. U.S.A. 98, 8950 (2001)]. In accordance with the experimental observations, we find that small 
nanoparticle concentrations induce an effective repulsion that prevents gelation caused by the in- 
trinsic van der Waals attraction between colloids. Higher nanoparticle concentrations induce an 
attractive potential which is, however, qualitatively different from the regular depletion attraction. 
We also show how colloid-nanoparticle size asymmetry and nanoparticle charge can be used to 
manipulate the effective interactions. 

PACS numbers: 82.70.Dd, 61.20.Ja 
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Colloidal particles find widespread application as pre- 
cursors for nanostructured materials, including advanced 
coatings, drug carriers, and colloidal crystals. These ma- 
terials are frequently fabricated from the liquid phase via 
hierarchical self-assembly processes that are governed by 
the interactions of the particles, their shape, and their 
size. The understanding and prediction of phase behav- 
ior and stability of a suspension relies on a fundamen- 
tal knowledge of the effective forces between colloids Q , 
which arise from a combination of direct interactions 
and indirect interactions mediated by the solvent and 
by other solute particles. Traditional methods to stabi- 
lize a suspension typically involve tuning of the effective 
interactions through charged groups or by grafting short 
polymer chains onto the colloidal surface. These mech- 
anisms, however, pose serious problems in certain sit- 
uations, such as the fabrication of close-packed colloidal 
crystals, where they cause an increase in the lattice spac- 
ing of the sedimented colloids and thus lead to cracking 
of the crystal upon drying. 

It is therefore of great fundamental and practical 
interest that an alternative stabilization strategy has 
emerg ed 0, Q. In these experiments, the addition of 
very small concentrations of charged zirconia nanoparti- 
cles was found to prevent the aggregation of near-neutral 
silica microspheres in aqueous suspension. Upon further 
increase of the zirconia concentration, microsphere aggre- 
gation was recovered, leading to a "window of stability." 
Lewis and co-workers Q ascribe the initial colloidal sta- 
bilization to the formation of a non-adsorbing "halo" of 
zirconia particles around the silica colloids, arising from 
the strong electrostatic repulsion between nanoparticles. 
However, while consistent with zeta-potential measure- 
ments that confirm a weak accumulation of nanoparticles 
near the colloidal surface, such halo formation is at vari- 
ance with the observation of stabilization at zirconia vol- 
ume fractions below 10 -3 , where the average nanoparti- 
cle separation is approximately 40nm, significantly larger 
than the electrostatic screening length k _1 « 2nm 



On the other hand, the reentrant gelation is attributed to 
the regular depletion attraction induced by the nanopar- 
ticles |4|. This is a potentially contentious argument as 
well, as it has been demonstrated that the depletion ef- 
fect can be substantially modified by non-additivity ef- 
fects 0, 0, which, in this case, would result from the 
electrostatic repulsions. 

Unfortunately, theoretical and computational ap- 
proaches experience severe difficulties in the treatment 
of this system. Analytical methods commonly incorpo- 
rate fluctuation and correlation effects only in an approx- 
imate manner. While simulations can explicitly account 
for these effects, inherent limitations in most computa- 
tional methodologies greatly restrict the range of acces- 
sible size ratios in multi-component mixtures. In this 
Letter, we overcome this problem by exploiting a novel, 
highly-efficient Monte Carlo scheme 8], which permits 
the explicit inclusion of interacting species of vastly dif- 
ferent sizes. Thus, we are able to determine the distribu- 
tion of nanoparticles and the colloidal potential of mean 
force as a function of nanoparticle concentration. Our 
findings constitute the first quantitative confirmation of 
the experimental observations and clarify the underlying 
physical mechanism of the stabilization. In addition, we 
provide explicit predictions for the role of two crucial pa- 
rameters, the magnitude of the nanoparticle charge and 
the colloid-nanoparticle size asymmetry. 

It is the purpose of our calculation to determine the ef- 
fective pair potential between colloidal particles induced 
by the nanoparticles. This effective interaction is inde- 
pendent of the intrinsic van der Waals attraction. If suffi- 
ciently repulsive it can counteract the latter and thus sta- 
bilize the suspension. Both species are modeled as hard 
spheres with diameter er m i cro and cr nan o , respectively. The 
aqueous environment (water and screening ions) is repre- 
sented as a homogeneous medium. In addition, we take 
into account the electrostatic double-layer interactions by 
means of the well-known Hogg-Healy-Fuerstenau (HHF) 
formula For two nanoparticles at a surface sepa- 
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ration D, this equation, under constant-potential condi- 
tions, reduces to 

►nano- nano 

*Lnoln[l + exp(- Kj D)] 
w £ e r 7rcr nano *^ ano exp(-/t£)) , (1) 

where the approximation is valid for kD ^> 1. $ nano 
is the zirconia zeta potential, eo the vacuum permit- 
tivity, and e r = 80 the dielectric constant of water. 
Since the silica spheres are negligibly charged, we ignore 
their electrostatic interaction. However, an electrostatic 
double layer interaction arises between the microspheres 
and the nanoparticles, which, for large size asymmetry 
a = cr micro /cr nano > 1 is described by [1] 
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where the approximation is again valid for kD 3> 1. 
While charge-regulating boundary conditions are even 
more appropriate than constant-potential boundary con- 
ditions for the interactions and J2Jl, both conditions 
have been found to lead to comparable results 01 • ^ n 
particular, the former can also yield an attraction be- 
tween neutral and charged surfaces . Indeed, whereas 
the presence of a microsphere-nanoparticle attraction is 
not considered explicitly in Ref . Q , both the supernatant 
measurements and the dynamic nature of the zeta po- 
tential measurement (which implies a cooperative move- 
ment of colloids and nanoparticles) indicate a certain de- 
gree of microsphere-nanoparticle association. It must be 
emphasized, however, that this attractive interaction by 
itself does not reduce the problem to a standard case 
of electrostatic stabilization resulting from an effective 
charge build-up on the colloids. Instead, as is shown be- 
low, the resulting "halo" is dynamic in nature and leads 
to effective interactions that are either attractive, repul- 
sive, or oscillatory, depending on nanoparticle concentra- 
tion. 

The presence of a short-range attraction is indeed con- 
firmed in Fig. ^ which provides a quantitative compari- 
son of the Monte Carlo simulations and the supernatant 
measurements. The adsorption of nanoparticles with di- 
ameter (T nano = 6nm around a microsphere with diameter 
o" micro = 0.6/im is determined as a function of nanoparti- 
cle concentration, at fixed volume fraction m i C ro = 0.1. 
In order to account for colloidal many-body effects and to 
handle the extreme size asymmetry (a = 100), we use the 
generalized geometric cluster algorithm Q for a system 
containing 50 colloids and up to 5 x 10 6 nanoparticles. 
In view of the small screening length and the presence 
of short-range hydration forces, we employ the approxi- 
mate expressions in Eqs. Q and (J2J) 14] . The amount 
of adsorbed nanoparticles, which form a layer with an 
approximate thickness <7 nano , is determined from the ra- 
dial distribution function integrated from contact to its 
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FIG. 1: Nanoparticle adsorption per colloidal microsphere 
as a function of zirconia nanoparticle concentration, at fixed 
^micro = 0.1. The open squares represent the simulation re- 
sults and the solid squares the experimental data 0. <^>nano = 
10 -3 corresponds to a concentration of 3650mg/l. 



first minimum. While the Monte Carlo data to some ex- 
tent overestimate the adsorption, both data sets exhibit 
a linear dependence on </> nano and a degree of adsorption 
that is far below 100%. The overestimation is compatible 
with the high zeta potential $ nano = 70mV, somewhat 
above the regime of validity of the Debye-Hiickel approx- 
imation employed in the HHF equations 0] . At volume 
fractions above 10~ 3 the mutual repulsion between the 
charged nanoparticles causes the adsorption to level off. 

Having established that the HHF equations model the 
experimental situation relatively accurately, we proceed 
to determine the effective (induced) colloidal interaction 
as a function of nanoparticle volume fraction. In the 
infinite dilution limit the potential of mean force fol- 
lows from the colloidal pair correlation function, V e s = 
— k^T In g{r). The simulations employ a cubic box with 
(f) micm — 0.002 16|. As shown in Fig. EJa), at very low 
0nano = 7 x 10 -4 a strong attraction exists at a separa- 
tion <7 nano , accompanied by a repulsive barrier of approx- 
imately 5fceT. The potential well arises from nanopar- 
ticles acting as "linkers" between two colloidal micro- 
spheres. However, upon increase of n ano to 10~ 3 the 
bridging attraction rapidly turns into a strong repulsion 
that extends over several nanoparticle diameters. Since 
the direct van der Waals interaction has a strength of only 
approximately — lksT at a surface separation er nano |3|], 
the total interaction exhibits a barrier. This repulsive 
barrier, which continues to increase with </> na no, is suffi- 
ciently strong to prevent gelation and thus leads to (ki- 
netic) stabilization of the suspension. The corresponding 
threshold agrees remarkably well with the experimentally 
observed gel-fluid transition near </> nano = 5x 10 -4 [2j. 

If the nanoparticle volume fraction is raised further, 
however, an attractive well arises for D w 4cr nano 
[Fig-EJb)] . As shown in the inset, the attraction strength 
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FIG. 2: Effective potential V^jk^T between a pair of col- 
loidal microspheres as a function of surface-to-surface separa- 
tion D. (a) At low nanoparticle volume fractions a bridging 
attraction is observed, which rapidly turns into a strong effec- 
tive repulsion, (b) At higher volume fractions an attractive 
minimum develops, with a strength that grows quadratically 

With </>nano (hlSet) . 



grows as 0n ano - At nano = 10~ 2 it reaches nearly 
— 3fceT, which is estimated to be sufficient to induce 
gelation. Again, this packing fraction agrees quite well 
with the experimentally observed reentrant gelation for 
0nano ~ 5 x 10~ 3 . The quadratic growth of the attrac- 
tion strength with </> n ano indicates that this is not a regu- 
lar depletion attraction, which has a linear concentration 
dependence Q. Indeed, examination of configurations 
shows only limited nanoparticle depletion between col- 
loidal pairs at D w 4<7 nano . 

In order to gain a microscopic understanding of the 
effective potentials we study the distribution of nanopar- 
ticles around a pair of colloids, for various separations of 
the colloidal particles. For clarity in the graphs, we em- 
ploy a smaller size ratio a — 40 (i.e., cr micro = 0.24/im) 
with </> n ano = 0.003, for which the effective interaction 
is strongly repulsive [cf. Fig. f4Ja) for nano = 0.001]. 
At a surface separation D > 5o- n ano [Fig. EUa)] the 
halos around both colloids are uniform, in agreement 
with a negligible effective interaction. If the colloids 



FIG. 3: [color online] Nanoparticle density distribution 
around two colloidal particles at a surface-to-surface dis- 
tance D. (a) At larger separations (D = 5.0cr nano ) the halos 
around each colloid are uniform, (b) Upon closer approach 
(D — 1.8<T nano ), where there is a strong colloidal repulsion, a 
depletion of nanoparticles in the gap region is observed. 



are forced closer together, the mutual repulsion between 
the nanoparticles leads to a redistribution of particles 
in the halo. At a separation D = 1.8cr nano [Fig. CJb)] 
most nanoparticles have been expelled from the gap re- 
gion and a depletion zone appears. While reminiscent 
of the depletion effect in hard-sphere mixtures the 
effective colloidal interaction in this case is strongly re- 
pulsive rather than attractive. This indicates that the en- 
ergy penalty resulting from nanoparticle rearrangements 
dominates over depletion-like entropic effects. 

The "haloing" stabilization mechanism clearly relies on 
rather generic features and hence may become of consid- 
erable practical importance. Since applications hinge on 
an understanding of the role of size and charge asymme- 
try, we have investigated a variety of parameter combina- 
tions. Figure flja) displays the repulsion resulting from 
a fixed volume fraction nano = 0.001 for colloids of di- 
ameter 0.60/im, 0.36/im, and 0.24/im (size ratio a = 100, 
60, and 40, respectively). Consistent with experimental 
observation Q , the effective colloidal repulsion increases 
with a. As shown in the inset, the dependence on size 
asymmetry is almost perfectly linear. The role of the 
nanoparticle charge is certainly complicated, as it affects 
both halo formation and the halo-halo interaction. In- 
deed, if the nanoparticle distribution around a pair of 
colloids at surface separation D is denoted by p(r';D), 
where r' is measured with respect to the first colloid, then 
the effective force on that colloid is 



F(D) = -jp(r';D)JLv mi 



micro — nano 



(r')dv' . (3) 



The nanoparticle charge will affect both p(r';D) and 
Vmicro-nano- In the absence of experimental data, we have 
performed simulations for the a — 100 case with a lower 
nanoparticle surface potential, $ nalM = 50mV [Fig.f4Jb)]. 
The effective interactions are found to be similar to those 
for higher surface potentials, with a comparable depen- 
dence on </> nano , but almost ten times larger volume frac- 
tions are needed to achieve the same effective colloidal 
repulsion. This emphasizes the important role of the 
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FIG. 4: (a) Effect of nanoparticle-colloid size asymmetry a on 
the effective colloidal interactions. At fixed nanoparticle size 
and volume fraction ^> n ano = 0.001, the repulsion increases 
with colloid size. As shown in the inset, the dependence is 
virtually perfectly linear, (b) Effect of nanoparticle charge. 
All parameters are identical to those in Fig. |2] except for the 
nanoparticle surface potential, which has been lowered from 
70mV to 50mV. As a result, much larger nanoparticle vol- 
ume fractions are required to achieve an appreciable colloidal 
repulsion. 

nanoparticle charge. Finally, we have performed control 
simulations in which the colloid-nanoparticle attraction 
is omitted. Neither the observed adsorption (Fig. nor 
the effective repulsions (Fig.|5J can be reproduced unless 
the nanoparticle concentrations are much higher than in 
the experiment or the nanoparticle repulsion has a much 
longer range, which in turn would be incompatible with 
the actual screening length « er nano /3 3]. 

In summary, the mechanism for colloidal stabilization 
through nanoparticle "haloing" |2( has been clarified by 
means of direct simulation. The geometric cluster al- 
gorithm, which provides an acceleration of many orders 
of magnitude |8|, has been indispensable for this: At a 
nanoparticle volume fraction (f> nano — 0.007 the calcula- 
tion involves 7 x 10 6 nanoparticles, requiring a simulation 
of O(10 4 ) hours. We have shown that an induced colloid 
nanoparticle attraction contributes to the effective repul- 



sion, but that it by no means reduces the situation to sim- 
ple steric or Coulombic stabilization. Haloing stabiliza- 
tion occurs over a window of nanoparticle concentrations, 
bracketed by two dissimilar gel phases, neither of which is 
a regular depletion gel. Our calculations indicate that the 
effective colloidal interactions also sensitively depend on 
nanoparticle charge and colloid-nanoparticle size asym- 
metry, so that the haloing mechanism indeed opens up a 
variety of opportunities to tailor materials properties. 
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Note: Upon completion of this work we became aware 
of a study by Karanikas and Louis 01 > wn0 use different 
techniques, but arrive at similar conclusions. 
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